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We study the Dirac quasiparticles in d-dimensional lattice systems of electrons in the presence of 
domain walls (d = 1), vortices (d = 2), or hedgehogs (d = 3) of superconducting and/or insulating, 
order parameters, which appear as mass terms in the Dirac equation. Such topological defects have 
been known to carry non-trivial quantum numbers such as charge and spin. Here we discuss their 
additional internal degree of freedom: irrespectively of the dimensionality of space and the nature of 
orders that support the defect, an extra mass-order-parameter is found to emerge in their core. Six 
linearly independent local orders, which close two mutually commuting three-dimensional Clifford 
algebras are proven to be in general possible. We show how the particle-hole symmetry restricts the 
defects to always carry the quantum numbers of a single effective isospin- 1/2, quite independently of 
the values of their electric charge or true spin. Examples of this new degree of freedom in graphene 
and on surfaces of topological insulators are discussed. 



I. INTRODUCTION 



Massless Dirac fermions arise often in the effective the- 
ories of condensed matter. The Fermi surface reduces to 
just two points in one dimension (d), and graphene, tt- 
nux state, and their generalizations provide realizations 
in the dimensions two and three. [l[ p-wave and d-wave 
superconductors represent yet another class of important 
examples. Q In the presence of some externally or dynam- 
ically generated order parameter (OP) Dirac fermions 
may become gapped, with the gap appearing as a mass- 
term in the Dirac Hamiltonian. Indeed, in any isotropic 
and time-reversal invariant Dirac system in d = 2, such 
as graphene, it is known that there are precisely 36 such 
linearly independent masses Q, which stand for various 
components of four insulating and four superconducting 
OPs [J, [5j . A pair of these masses may either commute 
or anticommute, with the latter ones being particularly 
interesting: they add in squares in the expression for the 
total gap, and this way cooperate in the reduction of 
the energy of the filled Fermi sea. The examples of such 
pairs in graphene include the two phase components of 
any superconducting (SC) order, or two spin components 
of the Neel OP, or of the quantum spin-Hall (QSH) topo- 
logical insulator. In principle, there exist numerous other 
and more exotic possibilities for such compatible OPs, in 
which one mass may even be insulating while the other 
is superconducting. [3j 

It is well-known that the domain walls and the vor- 
tex configurations in the masses of Dirac quasiparti- 
cles in d = 1 and d = 2 may carry electrical charge 
and/or spin. H-U The simplest way to understand this 
is on the example of insulating masses, when the elec- 
tric charge represents a good quantum number: for spin- 
1/2 fermions there are then two zero- modes of the Dirac 
equation, which when fully (half) occupied yield an in- 
tegrated vortex charge of one (zero) and spin zero (1/2). 
On the other hand, if the vortex configuration is in the 
SC order, for example, the particle number is not con- 
served, and the vortex charge is zero on average. Never- 
theless, the existence of zero-modes of the Dirac Hamil- 



tonian always leads to a non-trivial internal structure of 
the vortex, which, however, is independent of its elec- 
trical and rotational properties. Here we discuss such an 
additional quantum number that generally emerges inside 
these topological defects, and which may be understood 
as their "isospin" . In particular, we investigate how the 
definition, value, and the physical nature of the isospin 
depends on the type of defect and the dimension of space 
in which it is embedded. Special cases of the isospin 
components have already been identified in graphene: a) 
a vortex in the easy-plane Neel order was claimed to con- 
tain the hard-axis Neel component [lOj, and b) the s-wave 
SC vortex was argued to develop either a density-wave 
(DW) [ll| or the QSH state [l2j in its core. We generalize 
here these early examples, and in doing so, in particular, 
reconcile the seemingly conflicting results or refs. [11] 
and [HJ quoted above. 

We set up and discuss the most general formulation of 
this problem, and find that, whereas the OPs that consti- 
tute the isospins are dependent on detail, their number 
and algebra is universal and larger than previously found. 
It is independent of the nature of the OPs twisted into 
the defect, as long as their representative masses anticom- 
mute. It is also independent of the dimension of the sys- 
tem, if the defect is a topologically non-trivial mapping of 
the system's boundary onto the OP's target space. The 
structure we reveal is thus identical for a domain wall in 
a single mass in d = 1 , 0, 0| a vortex configuration in two 
anticommuting masses in d = 2.[8l-[l2| or a hedgeho g in 
three mutually anticommuting masses in d = 3.[ll. Il3l4l5l | 

What we find is six possible linearly independent OPs 
in the core of the defect, which close the Clifford algebra 
(7(3,0) <g> C(3,0),[lf| but with only the averages of the 
linear combinations of the orders belonging to the same 
factor algebra as finite. This implies that the defects in 
question always carry a single effective isospin degree of 
freedom, with the quantum number of 1/2. This is to be 
contrasted with the values of the real spin and the elec- 
trical charge of the defect, which both, and as we will see, 
even simultaneously, can average to zero. The mathemat- 
ical basis of this result lies in the somewhat arcane theory 
of real representations of Clifford algebras. 17j The ap- 
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pearance of two potentially non-zero isospins, in particu- 
lar, follows from the quternionic nature of the representa- 
tions relevant to lattice fermions of spin- 1/2. On the ex- 
perimental side, the complementary order in the defect's 
core should be visible as a local mass-gap in the scanning 
tunneling microscopy. The existence of the isospin could 
lead to new ways of manipulation of the defects by the ex- 
ternal probes that would couple to this particular degree 
of freedom. For example, the OP for the Kekule bond- 
DW order parameter couples to the pseudomagnetic field 
that derives from the wrinkles in the graphene sheet the 
same way the usual electric charge density couples to the 
real magnetic field. The possibility that defects in one 
OP carry with them quantum numbers of another is also 
suggestive of an exotic quantum phase transition from 
their proliferation. |18| 

A reconsideration of the defects in the insulating OPs 
becomes particularly interesting in light of our results, 
which imply that one of the potential isospins consists 
of a single insulating and two superconducting compo- 
nents. The average electric charge of such a defect is a 
continuous variable between zero and one. The chemical 
potential in this case polarizes the isospin in the direc- 
tion of the insulating component. It therefore causes the 
isospin to precess, which translates into a periodic time- 
dependence of the local superconducting phase. 

This paper is organized as follows. In the next section 
we define the Dirac-Nambu Hamiltonian in two dimen- 
sions, and discuss the consequences of the particle-hole 
symmetry, particularly of its antilinearity. The vortex 
Hamiltonian is then studied in Sec. III. In particular, in 
this section we show that the maximal Clifford algebra 
that contains the two gamma- and two mass-matrices 
appearing in the vortex Hamiltonian is closely related 
to C(2,5). A generalization to other dimensions and 
the connection to Wess-Zumino-Witten term is briefly 
described in Sec. IV. The way the isospin develops a 
non-zero expectation value from the zero-modes, and its 
further unitary evolution in time are considered in Sec. 
V. Some examples are provided in Sec. VI. The con- 
nection between the isospin on one side, and the usual 
spin and the electric charge on the other is discussed in 
the closing section, where we also summarize the main 
results. In the Appendices we provide the derivation of 
the Majorana representation used throughout the paper 
and give an elementary discussion of the crucial mathe- 
matical fact about the Clifford algebra C(2, 5), that is its 
quaternionic nature. 



have the standard form 



H = ctiki + a 2 fc 2 + 0(k 2 ), 



(1) 



where cti 



1,2 are eight- dimensional Hermitian ma- 



trices satisfying the anticommutation relation {a,, ctj} — 



25 a, and k 



-zV is the usual momentum operator. 



The dimension eight may be understood as arising, in 
graphene, for example, from two sublattices, two valleys, 
and two spin projections, but it is in fact required more 
generally by the time reversal invariance, as discussed in 
detail in ref. II 51 ] - To give a unified description of both 



the insulating and the superconducting mass-gaps one 
must double the number of fermionic components, in- 
clude both particles and holes into the Dirac spinor, and 
consider the sixteen- dimensional massive "Dirac-Nambu" 
Hamiltonian 

H = H k + mM (2) 
where, as usual, the kinetic energy part is 



H k = H ® {-H*) = Tih 



(3) 



and the matrices M and Ti, i = 1,2 are Hermitian. 
The symbol "0" stands for an orthogonal sum. To 
open a mass-gap we require that {M,T{\ = and 
M 2 = 1, so that the spectrum of H is the usual rela- 
tivistic ±\/k 2 + m 2 for m uniform. 

The Dirac-Nambu Hamiltonian H is by construction 
particle-hole symmetric: it anticommutes with the anti- 
linear operator I p h — (<ji ® I)K , where I is the eight- 
dimensional unit matrix, a are the Pauli matrices, and 
K is the complex conjugation. Since I 2 ^ = 1, a rep- 
resentation exists in which I p h — K [2(|, and H is 
imaginary. Explicit derivation is provided in Appendix 
1. Such a representation could be called "Majorana" , and 
we will assume it hereafter. In the Majorana representa- 
tion any mass-matrix M is also imaginary. Since, on the 
other hand, in the coordinate representation the momen- 
tum operator k is already imaginary, the two sixteen- 
dimensional matrices Ti , i = 1,2 in the kinetic energy 
are real. The separation of the matrices in the kinetic 
and mass terms into real and imaginary essentially de- 
termines the internal structure of the vortex, as will be 
discussed shortly. 



III. VORTEX AND ITS INTERNAL 
STRUCTURE 



II. DIRAC-NAMBU HAMILTONIAN 

To derive the isospin we take d — 2 and generalize to 
other spatial dimensions afterwards. Assume a space- 
and time-inversion invariant system of spin- 1/2 fermions 
on a lattice, with the Dirac point/points in its band struc- 
ture, such as graphene. The Dirac Hamiltonian for all 
low-energy excitations near the Dirac points will then 



Let us then assume a real-space vortex configuration 
in the mass term, 



M (x) = My = M l cos 9 + M 2 sin 9, 



(4) 



with (r, 9) as the polar coordinates, and m — m(r) in 
Eq. 2 and arbitrary. The only requirements are that 
the two Hermitian matrices Mi and Mi anticommutc, 
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TABLE I: The dimensions of irreducible real representations 
of Clifford algebras C(p,q), with p labeling the rows, and q 
labeling the columns. The normal font denotes the "normal" 
representations (with one Casimir operator), the emphasized 
font denotes the "almost complex" representations (with two 
Casimirs), and the bold font denotes the "quaternionic" rep- 
resentations (with four). 

{Mi,Mj} = 2Sij, again be imaginary, and anticommut- 
ing with the kinetic energy, i.e. {Mi,Tj} = 0. 

The spectrum of the Hamiltonian Hk + mMy now con- 
tains four zero-energy states. This follows from noticing 
that the matrices I\ and Mj, i = 1, 2 which appear in the 
Hamiltonian form a Clifford algebra C(4, 0), which, as 
is well known, has a unique four-dimensional irreducible 
representation. [2l| The Dirac-Nambu Hamiltonian can 
therefore be transformed into an orthogonal sum of four 
identical copies of the Jackiw-Rossi Hamiltonian, which 
contains one zero-mode in its spectrum. Q Any operator 
that anticommutes with the Hamiltonian may, depending 
on its occupation, pick up an expectation value precisely 
from those zero-energy states flOl - fia ] . To understand the 
internal structure of the vortex we therefore need to find 
the imaginary Hermitian matrices X that satisfy 

{X,T i }={X,M i } = 0, (5) 

i = 1,2, and X 2 = 1. 

In the Majorana representation the problem is solved 
by invoking some properties of the real representations of 
the Clifford algebras. To see how this algebraic structure 
emerges we first ask how many such mutually anticom- 
muting mass-matrices X exist. The answer is contained 
in Table I, where we extracted and summarized the re- 
quired information from ref . (l7| , both on the dimension- 
ality and the type of the irreducible real representations 
of the Clifford algebras C(p,q). Since 1^, iMi, and iX 
are all real, and Tf = +1 and {iMi) 2 = (iX) 2 = -1, 
we need to determine the maximal value of q so that 
for p > 2 the dimensionality of the real representation 
of C{p, q) is sixteen. An examination of the table im- 
mediately implies that the maximal q — 5, as observed 
already in the ref. Q. There are thus three additional 
mutually anticommuting masses: X = M3, M4, M5, so 
that the set (Ti,iMj), i = 1,2, j = 1,2,3,4,5 provides 
a real irreducible sixteen-dimensional representation of 
C(2,S). 




m; 

FIG. 1: The three mutually anticommuting OPs (the first 
isospin's components) in the defect's core. The OP at each 
corner commutes with the symmetry operator at the opposite 
side of the triangle, which rotates the remaining two OPs into 
each other. There also exists the second triangle of OPs, with 
Mi Mi (the second isospin's components) and Ei — > E[, 
which is not shown. All the matrices are imaginary, and all 
the matrices of the first triangle commute with those of the 
second. 

This is not all, however. The essential feature of 
the real representation of C(2,5) is that it is quater- 
nionic: besides the usual unit matrix, there exist three 
non-trivial linearly independent sixteen dimensional ma- 
trices which commute with the whole representation. 
This facilitates the construction of three additional solu- 
tions: X = M3, M4, M5. While mutually anticommuting 
{{Ml, Mj} = 26ij), these solutions commmute with the 
first three: [Mi, Mj] = 0, for i,j = 3,4,5. The mathe- 
matical details of this construction are provided in Ap- 
pendix 2. 

All possible "five-tuplets" of anticommuting masses 
which respect the particle-hole symmetry have been 
listed with the help of a computer in ref. Q. From 
this list it is evident, under close inspection, that each 
pair of symbols, irrespective of what physical OPs they 
may represent, can be found in exactly two different five- 
tuplets. Any set of three, or four, symbols, on the other 
hand, belongs to a unique five-tuplet. All of these facts 
follow from the representation theory as explained above. 



IV. OTHER DIMENSIONS AND 
WESS-ZUMINO-WITTEN TERM 

The six masses that we found to anticommute with 
the Dirac-Nambu Hamiltonian in d = 2 with a mass- 
vortex close the algebra C(3, 0) (g> C(3, 0), with the four- 
dimensional Hilbert space of zero-energy states Hq as its 
invariant subspace. Table I implies, however, that this 
result is independent of the spatial dimensions d: since 
the dimensions of the matrices appearing in the Dirac- 
Nambu Hamiltonian are 2 , d real T- matrices and 
d + 3 imaginary masses, upon multiplication of the latter 
by i, form a real representation of C{d, d + 3), which is 
always quaternionic. Assuming the topological defect in 
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d of the anticommuting masses again leaves the algebra 
of potential core order parameters C(3, 0) <g> C(3, 0). The 
general result can then be summarized as in Fig. 1. 

In any dimension d, the number of mutually anticom- 
muting masses, rf + 3, is precisely right for the emergence 
of the Wess-Zumino-Witten term upon the integration 
over fermions. [ill [HHUI The Berry phase of the three 
OPs that can form in the core of the defect of the re- 
maining d masses will thus match the one of a spin- 1/2- 
like object, as discussed in [24j], and in accord with our 
assignment of the isospin's quantum number. 



V. TWO ISOSPINS AND ZERO-MODE 
SUBSPACE 

A. State of isospin 

Since the zero-mode subspace Hq is four-dimensional, 
we may choose a basis in which the first three core OPs 
(Mi) are represented by Si, = Oj ® do, and the second 
three (M[) by S 2 i = <?a ® a% (i = 3,4,5). It would then 
seem that the averages of both, or possibly none, of the 
two isospins could be finite in the ground state. That this 
however, is impossible, is guaranteed by the antilinearity 
of the particle- hole symmetry in a simple but rather inter- 
esting way. Consider a perturbed Dirac-Nambu Hamil- 
tonian 

H = H k + mM v + eO, (6) 

where is an operator so far neglected, which derives, 
for example, from the underlying lattice, magnetic field, 
finite chemical potential, or a similar perturbation. For 
small enough parameter e we may neglect the mixing 
of the zero-energy and finite-energy modes, and to the 
leading order consider only the resulting split of the zero- 
modes. We now show that the result of any such splitting 
can only be that either (Si) ^ and (S2) = or vice 
versa. 

Since in the Majorana representation the matrix O for 
any perturbation is also imaginary, O anticommutes with 
some antilinear operator A that represents the particle- 
hole symmetry within the zero-mode subspace Hq. In- 
troduce the basis in Ho: {|±,±)}, with S a 3\si,s 2 ) = 
s a |si, S2), where s a — ±, and a — 1, 2 labels the isospin. 
If is an eigenvector of the matrix representing O in 
"Ho, its eigenvalue has the opposite sign of the eigenvalue 
of the eigenstate A\if>). The two states are thus orthogo- 
nal. Since S a , a = 1,2 are also odd under particle- hole 
symmetry, it follows that within Hq A is formally equiv- 
alent to the time reversal of spin. In general then, if 

M = a\+,+)+b\+,-) + c\-,+)+d\-, ->, (7) 

its particle-hole partner is 



where A|±) = ±i|=p). The orthogonality of the two states 
implies that 



which is precisely what disentangles \ip) into a product 
state: 

W = (l+)i + -|->i)(«l+>2 + 6|-) 2 ) = \ft)i\m) 2 , (10) 
a 

where, in the usual notation, 

**• § a $ a = |tJ Q (11) 

for a — 1, 2, and t 2 = 1. 

If the state \ip) = \n)i\m)2 corresponds to a negative 
eigenvalue of O, the state A\ip) = \ — f?)i| — rh)^ cor- 
responds to a positive eigenvalue. Since the other two 
eigenstates are also the product states, and need to be 
orthogonal to the first two, there are just two options: it 
is either that the state | — n) 1 \m) 2 has then the negative 
energy, and \n)i \ — 771)2 positive, or vice versa. In the 
former case, (S2) — m, and (Si) = 0. The state of the 
isospins may thus be understood as a mixture, given by 
the statistical operator (i.e. the density matrix) 

P = iffo® \m)(fn\. (12) 

In the latter case, of course, = n, and (S 2 ) = 0. 

The state of the two isospins is a product of the maximal 
mixture (1/2<to) for one, and the pure state (|m)(m|) for 
the other isospin. 

B. Time evolution 

The unitary evolution cannot exchange the two 
isospins: if (Si) = t = 0, for example, it remains zero 
at t > 0. To see this, assume that the isospin's state ini- 
tially is p(0) = (l/2)i7o <8> p 2 (0) with a general statistical 
operator p 2 (0). It evolves according to the Schrodinger 
equation 

ip(t) = [0,p(t)]. (13) 

The particle-hole symmetry implies that in general the 
perturbation, which may be time-dependent, within the 
zero-mode subspace has the form 

O = f(t)Pi(t) ®cr + h(t)a ® P 2 (t), (14) 

where Pj (t) = \ fi(t)){ri(t)\ are time-dependent projectors 
on some one-dimensional subspaces of the first and the 
second isospin, and f(t) and h(t) are arbitrary functions 
of time. The time-dependent state is therefore p(t) — 
(l/2)oo ® Pi{t), with 



A\iP) = -a*\-,-) + b*\-,+)+c*\+,-)-d*\+,+), (8) 



ip2(t) = h(t)[P 2 (t),p 2 (t)], 



(15) 
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and, as evident, (Si) — at all times. 

The topological defects thus ultimately carry with 
them a single non-trivial isospin-1/2, of the physical 
size of the spatial extent of the zero-modes, which is 
~ l/m(oo). The OPs forming the components of the 
isospin depend on the system. Let us consider some ex- 
amples next, starting from the most familiar. 

VI. EXAMPLES 

A. Superconducting vortices 

Let us consider graphene (d — 2), as an example. If 
the vortex configuration is in the two components of the 
s-wave SC, core's isospin components are the charge-DW 
and the two singlet bond-DWs,[ll| or the three compo- 
nents of the QSH state The real spin in the former 
case is zero. If the superconductor is a spin singlet, as 
in this example, one of the isospins must be a vector 
under rotations, since all three spin-rotation generators 
commute with the Hamiltonian, and will appear as the 
Casimir operators £3, £4,^5 in Fig. 1, for one of the 
isospins. Obviously, in this case it is the QSH state that 
carries the (true) spin 1/2, which is why it is favored by 
the Zeeman term in presence of the magnetic field. [12] In 
general, if the matrices Mi and Mi in Eq. (4) represent 
two phase components of some superconducting order, [5j 
all six core's OPs are insulating. The electric charge is 
not a good quantum number, but on average it always 
vanishes, independently of the state of isospin. The lat- 
tice itself, for example, favors the bond-DW in the core: 
this is thus an example of a vortex with vanishing both 
(average) charge and true spin, and with only the isospin 
left as finite. 

The situation changes if one twists two spin compo- 
nents of a triplet SC; x and y components of the f-wave 
SC in graphene, [!, [25| for example. Each potential 
isospin then has a single insulating and two supercon- 
ducting components: 1) CDW, s-wave SC, z component 
of the f-wave SC, and 2) z components of Neel, and two 
Kekule triplet superconducting states. In either case the 
vortex core may stay superconducting by simply chang- 
ing the type of the superconducting state, i. e. by form- 
ing a texture. Neither electric charge nor spin are good 
quantum numbers in this case, but the isospin is well- 
defined. 



B. Insulating vortex 

Next, consider a vortex configuration in two insulating 
OPs, when the particle number, and hence the electric 
charge, is a good quantum number. For instance, the 
spin-singlet (Kekule) bond-DW @ is favored by the cou- 
pling to in-plane optical phonon at the Dirac point in 
graphene, and has been prop osed as a candidate for the 
OP in the magnetic field. |26l428j | The isospin components 



in the vortex configuration in this case are, 1) the three 
components of the Neel OP, or 2) the charge-DW and 
the two phase components of the s-wave SC. The vortex 
is (spin) rotational symmetric, and the spin is 1/2 and 
zero, respectively. The electric charge for the first isospin 
is zero. 

Before discussing the electric charge of the second 
isospin, let us consider the vortex formed out of the x 
and y components of the Neel OP, favored by the Hub- 
bard on-site repulsion. [29| One set of isospin components 
contains the z-component of the same Neel order [l(| 
and two spin-singlets bond-DWs, and the other consists 
of the z-component of the QSH state and the two phase 
components of the z-component of the f-wave triplet SC 
Spin is not a good quantum number. 

In both examples introduced above, one isospin (call 
it Si) is purely insulating, whereas the second one (S2) 
has one insulating and two superconducting components. 
It is easy to see that this is in general true; the particle 
number is conserved, and thus the number operator must 
be one of the Casimir operators, let's say E3 in Fig. 1. 
The masses M' 4 and M' b must then stand for two phase 
components of some superconducting order, rotated into 
each other by a transformation generated by the number 
operator. 

Let us further assume that in the presence of 
some perturbation we found (S2) — n, where n — 
(sin 9 cos cj>, sin 9 sin <j>, cos 9) , with S23 — Co ® °3 as the 
insulating component of the second isospin. The particle 
number operator within Jig is then also N — <7o <£> 03, 
and the average of the vortex electric charge is (N) — 
(fi\as\n) = cos9. Of course, in this example the charge 
will be sharp only if 9 = 0, which will be the case if 
the perturbation in question is the chemical potential 
/j,N, for instance. In the ground state then 9 — 0, and 
the charge is back to unity. One can imagine, however, 
that a proximity to another superconductor could "po- 
larize" the core so that 9^0. At a finite chemical po- 
tential the superconducting phase will then oscillate as 
4>(t) — 0(0) + /-it, whereas 9 is constant, just as in the 
precession of an ordinary spin. 



C. Topological insulator 

Our results do not, however, readily apply to a super- 
conducting vortex on an isolated surface of a topological 
insulator, if the original fermion has only two compo- 
nents. After the particle-hole doubling the dimension of 
the matrices involved becomes four, and Table I implies 
that for p — 2 the maximal value of q is only two. The set 
of the core's order parameters is simply empty, in accord 
with the single zero-energy Majorana state being pinned 
to the chemical potential. [3(| If we consider a thin film 
geometry with two possibly coupled surfaces, however, 
the dimension of the Dirac-Nambu Hamiltonian doubles 
to eight, and the maximal value of q increases to q = 3. 
The situation still differs qualitatively from the one we 
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studied hitherto in that there is now a single isospin in 
the core, and with the discrete Ising, and not Heisenberg 
symmetry. If the vortex is in the s-wave SC, for example, 
one finds that the core's Ising isospin is the time-reversal- 
invariant exciton of ref . 13 ill . 



VII. SUMMARY AND DISCUSSION 

To conclude, we showed that some typical topological 
defects in the presence of Dirac fermions, besides their 
usual electric charge and true spin, which, depending on 
the defect, may or may not be finite, always possess a 
non-zero isospin. The components of the isospin are the 
local masses, i. e. order parameters, that develop a fi- 
nite expectation value in the center of the defects. While 
in principle two different such isospins are always possi- 
ble, only one can have a non-trivial quantum number in a 
given physical situation. This new degree of freedom cou- 
ples to different external probes, and in principle could 
be manipulated similarly to the electron spin. 

The connection between the isospin, discussed here, 
and the electrical charge and the true spin of the vortex 
now becomes transparent. A particular split of the zero- 
modes favored by the perturbations to the Dirac Hamil- 
tonian will result in a component of one of the two pos- 
sible isospins developing a non-zero average value; say, 
for instance, (wM^) ^ 0. The vortex is thus by this 
mechanism spontaneously turned into a meron. There is 
a unique non-trivial symmetry generator that commutes 
both with the vortex Hamiltonian, and with the matrix 
M3: it is E3 in Fig. 1. The commutation relations be- 
tween the six masses M, and M[ and the six Casimir 
operators Ei and E[ imply, however, that within the zero- 
mode subspace, E3 and M' z are represented by the same 
matrix, and therefore it is also (^^E^) ^ 0. The ma- 
trix E3 can be the number operator, in which case the 
latter average would be proportional to the density of the 
electric charge. It can also be the generator of rotations 
of the electron spin, in which case the vortex will have 
the spin 1/2. The examples of both were discussed in 
the previous section. But E3 can also be neither. For 
a singlet-SC vortex on graphene's honeycomb lattice, for 
example, if the non-zero isospin is Kekule bond-DW, E3 
will be one of the generators of valley rotations. 
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IX. APPENDIX 1: MAJORANA 
REPRESENTATION 

Particle-hole symmetry dictates that the mass matrix 
M satisfies the constraint: 

M = -(a 1 ®I)M T {a 1 ®I) (16) 

where I is the eight-dimensional unit matrix, and a are 
the Pauli matrices. The constraint is the expression 
of the Fermi statistics, and it implies that the bilinear 
* t M* is finite. 

The determination of all possible masses is therefore 
made non-trivial by the constraint. Luckily, it can be 
simplified by a unitary transformation M = UMU\ so 
that M would satisfy M = -M T . Here U = U 2 <E)I, with 
the two-dimensional unitary matrix U% satisfying 

tfo = U 2 a 1 U 2 T . (17) 

The general solution is then 

U 2 = (iO^V^V^V^- 1 ^ 3 . (18) 

4> = 7r/4, for example, yields the matrix of ref. [20| . 
Since the matrix M is also Hermitian, it follows that 
M = —M*, i. e. that it is purely imaginary. Also, 
{M, H k } = 0, where H k = UH k W. 

Writing on = Re(a)i + ilm(a)i, the kinetic energy of 
the Dirac-Nambu Hamiltonian takes the form 

H k = (fo <8> Re(a)i + ia 3 (gi Im(a)i)ki = Tih. (19) 

Since, on the other hand, U 2 a^lj\ = o~ 2 , after the unitary 
transformation the same kinetic energy becomes 

Hk = (o- ® Re{a)i + ia 2 ® Im{a)i)k. t = tiki, (20) 

with the matrices Ti being manifestly real. 

X. APPENDIX 2: QUATERNIONIC 
REPRESENTATION OF C(2,5) 

The crucial feature of the real sixteen-dimensional rep- 
resentation of the Clifford algebra C(2,5) is that it is 
quaternionic: besides the unit matrix, there exist three 
additional real Casimir matrices, Ki, i = 3,4,5, which 
obey the quaternion algebra, 

KiKj = -Sij + e ijk K k , (21) 

and which commute with all seven matrices {Ti,iMj}, 
i = 1,2, j = 1,2, 3, 4, 5. Defining then the imaginary ma- 
trices as Ei = iKi, i = 3, 4, 5, we see that they mutually 
anticommute, {Ei,Ej} = 2Sij, and close the SU(2) Lie 
algebra. Their existence enables one to construct addi- 
tional solutions for the matrices X as follows. Consider, 

M[ = E l T 1 T 2 M 1 M 2 (22) 
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and i = 3,4, 5. Obviously, {M-, Vj} = {M[, Mj} = for 
j = 1,2. Mj, i = 3,4,5 are both imaginary and Hcrmi- 
tian, and {M[, Mj} = 2<%. The set {Mi, Mj}, i = 1,2, 
j = 3,4,5 is then a five-tuplet of mutually anticommut- 
ing masses, different from {.M,Mj}. We also see that 
[Mi, Mj] = for i, j = 3, 4, 5, and, by Eq. 1, 

iE 3 Mi = M5, (23) 

together with cyclic permutations of the indices. M/, 
i = 3, 4, 5 therefore form a vector under transformations 
generated by Ei. Of course, the real representation of 
C(2, 5) provided by {T u iM 3 ,iM' k }, i, j = 1, 2, fc = 3, 4, 5, 
is also quaternionic, with its own set of Casimir operators 
E[, i = 3, 4, 5, so that it is also true that 

M l = E' l Y 1 Y 2 M 1 M 2 . (24) 

The origin of the non-trivial Casimir operators is easy 
to understand. The maximal number of complex sixteen- 
dimensional mutually anticommuting matrices is nine. 
If they are all to square to +1, five of them, call them 
Ri, i = 1,2,3,4,5 can be real, and the other four, I i} 
i = 1,2,3,4, imaginary. This is a direct generalization 
from the case of two-dimensional Pauli matrices, and is in 
accord with the normal nature of the real representations 
of C(q + l,q) in the Table I in the text. We may take 



that, for example, 

4 

R 5 = Y[RiIi. (25) 

Form then the subset of six matrices {Ri, Ik}, i = 1,2, 
and k = 1, 2, 3, 4. Obviously, we can construct an unique 
additional imaginary matrix out of the three omitted real 
matrices as 

h=iR 3 R 4 R 5 , (26) 

which then anticommutes with all six matrices of the 
subset. The three imaginary Casimir operators are then 
simply 

Ei = 2 e i3kRjRk, (27) 

where i, j, k — 3, 4, 5 in the last equation. Generalization 
to other dimensions is straightforward. One can similarly 
understand the "almost complex" real representations of 
C(p,p+1), with two Casimir operators, which is relevant 
for the vortices on two surfaces of topological insulator 
thin films discussed in the text, or for spinless fermions 
on honeycomb lattice. 
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